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Abstract 

Coronaviruses from both the Alpha and Betacoronavirus genera, interfere with the type I interferon 
(IFN) response in various ways, ensuring limited activation of the IFN response in most cell types. 
Of Gammacoronaviruses that mainly infect birds, little is known about activation of the host 
immune response. We show that the prototypical Gammacoronavirus, infectious bronchitis virus 
(IBV), induces a delayed activation of the IFN response in primary renal cells, tracheal epithelial 
cells and in a chicken cell line. IfnB expression in fact, is delayed with respect to the peak of viral 
replication and accompanying accumulation of dsRNA. In addition, we demonstrate that MDAS is 
the primary sensor for Gammacoronavirus infections in chicken cells. Furthermore, we provide 
evidence that accessory proteins 3a and 3b of IBV modulate the IFN response at the transcriptional 
and translational level. Finally, we show that, despite the lack of activation of the IFN response 
during the early phase of IBV infection, signalling of non-self dsRNA through both MDAS and TLR3 
remains intact in IBV-infected cells. Taken together, this study provides the first comprehensive 


analysis of host-virus interactions of a Gammacoronavirus with avian innate immune responses. 


Importance: Our results demonstrate that IBV has evolved multiple strategies to avoid activation 
of the type I interferon response. Taken together, the present study closes a gap in the 
understanding of host-IBV interaction, and paves the way for further characterization of the 
mechanisms underlying immune evasion strategies as well as pathogenesis of 


Gammacoronaviruses. 
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Introduction 


Coronaviruses constitute a large family of positive-stranded RNA viruses and cause a range of 
human and veterinary diseases. Infectious bronchitis virus (IBV) is the prototype avian coronavirus 
from the Gammacoronavirus genus and the causative agent of a highly contagious respiratory 
disease of major economic importance to the poultry industry (1). IBV enters the avian host 
through the respiratory tract, where it causes destruction of the epithelium leading to respiratory 
distress and initiation of secondary bacterial infections. Depending on the strain, IBV can also 
spread to other epithelial surfaces such as the gastrointestinal tract, the kidneys and the oviduct, 
the latter causing problems in egg production and quality (1-6). Contrary to coronaviruses from the 
Alpha and Beta genera, including human coronavirus HCoV-229E, Severe Acute Respiratory 
Syndrome (SARS-CoV), Middle East Respiratory Syndrome (MERS-CoV) and mouse hepatitis virus 
(MHV), very little is known about how Gammacoronaviruses including IBV evade or interfere with 


innate immune responses of their host. 


Innate immune responses consist of a network of antimicrobial mechanisms, of which the type 1 
interferon (IFN) response is an essential defence mechanism against viruses. Typically, the type ፤ 
IFN response, from hereafter referred to as IFN response, is initiated upon activation of host 
pattern recognition receptors (PRRs), present in all animal cells. Two families of PRRs have been 
shown to be involved in the recognition of RNA viruses namely the membrane-bound Toll-like 
receptors (TLRs) and the cytosolic RIG-I-like receptors (RLRs) (7). The primary ligands for 
activation of these PRRs are double-stranded RNA (dsRNA) and 5' triphosphate-containing RNA, 
normally absent from uninfected host-cells. Activation of RLRs leads to the transcription of genes 
encoding type I interferons (IFNa and IFNB). These interferons are secreted from the infected cell 
providing a signal for the infected as well as the neighbouring cells that induce the transcription of 


anti-viral effector genes collectively called interferon stimulated genes (ISGs). 
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The ability of a virus to replicate and produce infectious progeny depends for a large part on its 
ability to avoid induction or counteract the IFN response of its host. Indeed, a common feature of 
Alpha- and Betacoronaviruses, including HCoV-229E, SARS-CoV, and MHV, is their limited 
activation of the IFN response (8-13). This limited activation can be partially explained by 
intracellular membrane rearrangements that might shield dsRNA and other viral components from 
recognition by host PRRs (14, 15). In addition, coronavirus nsp16 displays 2'-O-methylase activity, 
which results in 2'-O-methylation of a ribose moiety on the 5’ cap of coronavirus mRNAs, making 
them indistinguishable from host mRNAs (16). Furthermore, many other coronavirus proteins, such 
as nspl, nsp3, the nucleocapsid and many of the accessory proteins have been shown to interfere 


with the IFN response in various ways (reviewed in (17, 18)). 


Interaction between Gammacoronaviruses and innate immune responses of their avian hosts is 
poorly understood. Early studies on Gammacoronaviruses in chicken suggest that IBV-induced IFN 
production is variable and dependent on both virus strain and cell type. (19-22). Further, two 
transcriptional studies on tissues collected after in vivo and in ovo IBV infections, found only limited 
upregulation of ISGs at 1 - 3 days post-infection (23-25). Functional studies using IBV Beaudette 
showed that it induced cell-cycle arrest and apoptosis (26, 27), that IBV interacts with elF3f (28) 
and that IBV inhibits protein kinase R activation, thereby maintaining protein synthesis (29). 
Although these studies did provide a number of details on the interactions between IBV and the 
host cell, most experiments were carried out in Vero cells. This non-avian cell line is one of the 
very few cell lines in which the IBV-Beaudette strain has been adapted to grow, facilitating in vitro 
experiments. Vero cells, however, lack the IfnB gene, preventing them from mounting a type I IFN 
response (30, 31), reducing the value of Vero cells for research on innate immune responses to 
IBV. In addition, the Beaudette strain is non-pathogenic in vivo with limited replication in host 
tissues (32), reducing the value of these in vitro studies for translation to in vivo situations. For 
these reasons, we used pathogenic isolates of IBV to infect primary chicken cells, and a chicken cell 


line, as these isolates are known to infect, spread and cause clinical disease in vivo. 


In the current study, we show that IBV infection leads to a significant induction of IfnB transcription 
through an MDA5-dependent activation of the IFN response, albeit delayed with respect to both 
virus replication and accumulation of dsRNA. This delayed induction of IfnB was further confirmed 


through RNA FISH analysis showing that accumulation of Ifn8 mRNA is restricted to IBV-infected 
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and not neighbouring uninfected cells. Although the time lag between accumulation of dsRNA and 
induction of IfnB transcription might suggest that IBV interferes with recognition of dsRNA, we 
observed that sensing of exogenous (non-self) dsRNA remained functional in IBV-infected cells. 
Using mutant IBV viruses we demonstrate that both accessory proteins 3a and 3b are involved in 
limiting IfnB expression, as both 3a and 3b null viruses induced increased IfnB expression. 
Nevertheless, 3a and 3b seem to have a differential effect on IFN protein production, infection with 
3a null virus induced lower IFN levels whereas a 35 null virus increased IFN production compared 
to the parental virus. Altogether, our data suggest that IBV delays but does not prevent detection 
by MDA5, and that accessory proteins 3a and 3b modulate the IFN response in avian cells. This is 
the first study addressing immune evasion and interference strategies of IBV in chicken and not in 
mammalian cells, providing information essential to further understanding of the pathogenesis of 


Gammacoronaviruses. 


MATERIALS AND METHODS 


Cells 

Chicken embryonic kidneys were aseptically removed from 17- to 19-day-old chicken embryo's 
(Charles River, SPAFAS). A cell suspension was obtained by trypsinisation for 30 min at 37°C and 
filtered through a 100 uim mesh. The resulting chicken embryo kidney (CEK) cells were seeded at 4 
x 10? cells/cm? in 199 medium (Invitrogen) supplemented with 0.596 fetal bovine serum (FBS, 
FASC) and 1% PenStrep (Gibco®, Invitrogen). Chicken trachea cells were isolated from 8- to 10- 
week-old chickens (white leghorn). Tracheas were collected in ice-cold PBS, washed and stripped 
from adipose tissue. Trachea were filled with a solution of 3.5 U/ml protease type XIV (Sigma), 4 
U/ml DNase I (Qiagen) and 1% PenStrep in EMEM, sealed with clamps, and incubated overnight at 
49C. The next day, cells lining the luminal side of the trachea were flushed out with cold EMEM, 
filtered through a cell strainer and seeded at 4 x 10? cells/cm? in DMEM supplemented with 10% 
FBS and 1% PenStrep. The RIG-I"*, RIG-I“°, MDA5"' and MDAS'? MEFs were provided by Prof. S. 
Akira (33). The MAVS"' and MAVS'? MEFs were provided by Z.J. Chen (34). DF-1, CEC-32 and 
MEFs cells were cultured in DMEM (Gibco®, Invitrogen) supplemented with 10% FBS and 1% 


PenStrep. All cells were incubated in a humidified incubator at 37°C and 5% CO;. 
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DF-1 IfnB-luc reporter cell line 

DF-1 cells were transfected using Fugene (Promega) according to the manufacturer's instructions 
with a construct expressing firefly luciferase under the control of the -110bp proximal region of the 
human IFNB promotor (35). Stably expressing cells were selected over a period of 3 weeks using 
geneticin (500 pg/ml). DF-1 IfnB-luc stable cells were cultured in DMEM supplemented with 10% 


FBS and 1% PenStrep and were not further subcloned. 


Viruses 

IBV-M41, IBV-QX and IBV-Italy-O2, Rift Valley Fever Virus clone 13 (RVFV Cl13) and Infectious 
Pancreatic Necrosis Virus (IPNV) were obtained from Merck Animal Health, Boxmeer, The 
Netherlands. Sindbis-GFP was a kind gift from J. Fros, (Laboratory of Virology, Wageningen 
University). IBV Beaudette, strain Beau-R, as well as the generation of the ScAUG3a, ScAUG3b, 
ScAUG3ab and ScAUG5ab Beau-R null viruses has been published previously (36-38). In these 
mutant IBV viruses, the start codons of the indicated accessory genes were mutated to stop 
codons. All IBV strains were amplified and titrated on CEK cells. Sindbis-GFP was amplified on baby 
hamster kidney (BHK) cells and titrated on CEK cells. RVFV Cl13 was amplified and titrated on Vero 
cells, an African green monkey cell line. IPNV was amplified and titrated on the CHSE-214, 
Chinook-salmon cell line. IPNV was inactivated by 20 min UV exposure on a 48W BXT-26-M 


(Uvitec). 


Poly I:C stimulation and RNase treatment 

Polyinosinic-poly(C) [p(I:C)] sodium was purchased from Sigma, dissolved in nuclease-free water 
and stored at -80°C. p(I:C) was either directly added to the medium or transfected using 
Lipofectamine 2000 (Invitrogen) as per manufacturer's instructions. DF-1 cells (3 x 10"/well) were 
cultured in 24 well plates and transfected with 500 ng p(I:C). RNase treatment of CEK cell culture 
supernatant was performed by addition of 10 pg/ml RNase A (Invitrogen) before IBV infection or 


before stimulation with 2 ug/ml p(1:C). 


RNA isolation and cDNA synthesis 
Approximately 8 x 10" CEK cells or 3 x 10° DF-1 cells were lysed in RLT buffer (Qiagen) at various 
time points after treatment or infection. RLT cell lysis buffer was spiked with 1 ng/sample of 


luciferase mRNA (Promega) immediately prior to RNA isolation. Luciferase expression will later be 
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used as external reference gene for normalization during the gene expression analysis. Total RNA 
was isolated using the RNeasy Mini Kit (Qiagen) according to the manufacturer's instructions, 
including an on-column DNase treatment with RNase-free DNase (Qiagen). Before cDNA synthesis 
of 0.5-1 ug total RNA, a second DNase treatment was performed using DNase I, amplification 
grade (Invitrogen). Synthesis of cDNA was performed using SuperScript III (Invitrogen) using 
random primers. cDNA samples were further diluted 1:50 in nuclease-free water before real-time 


quantitative PCR analysis. 


Gene expression analysis 

Real-time quantitative PCR was performed on a Rotor-Gene 6000 (Corbett Research), using 
Brilliant SYBR Green quantitative PCR (Stratagene) and primers (39-42) as listed in Table 1. Cycle 
thresholds and amplification efficiencies were calculated by the Rotor-Gene software (version 1.7). 
The relative expression ratio of the target gene was calculated using the average reaction efficiency 
for each primer set and the cycle threshold (C:) deviation of sample vs. control at time point Oh, as 
described in (43). For calculation of the fold change of IBV total RNA, C; deviation was calculated 
versus C, 30, as no IBV was present in the non-infected cells that were used as control in all the 
experiments. Because expression of various housekeeping genes was unstable during virus 
infections at time points later than 24 hours (data not shown), gene-expression ratios were 


normalised using an external reference gene (luciferase). 


Immunohistochemistry 

CEK cells were seeded on fibronectin-coated glass Biocoat coverslips (BD Biosciences) at a density 
of 1 x 10° cells/cm?. After incubation at 37°C for 48 hours, cells were infected with IBV strain M41 
at an MOI of 1, and fixed at different time points with 3.7% paraformaldehyde and permeabilized 
using 0.1% Triton X-100. Infected cells were probed with anti-dsRNA antibody (English & Scientific 
Consulting) and polyclonal chicken serum raised against IBV M41 was obtained from Merck AH. 
Detection was done performed using Alexa 488 goat anti-mouse antibody (Invitrogen) and FITC 
labelled goat anti-chicken antibody (Kirkegaard and Perry laboratories). Nuclei were stained with 
4',6-diamidino-2-phenylindole (DAPI). Cells were imaged using a Zeiss Primo Vert microscope and 


Axiovision software. Image overlays were made in Image]. 
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RNA fluorescence in situ hybridization 

RNA fluorescence in situ hybridization (FISH) was performed according to previously described 
protocols (44-46). ል set of forty RNA FISH probes (20 bp), each labelled with one CAL Fluor® Red 
610 fluorophore and targeting chicken IfnB (ENSGALT39477), was designed using the Stellaris® 
probe designer (Biosearch Technologies; https://www.biosearchtech.com/stellarisdesigner/). The 
coding sequence of chicken Ifnf is 601 bp, therefore to accommodate the optimum number of 
fluorescent probes (48, as explained in reference (44)), the 3'UTR was included in the probe design 
tool. CEK cells were grown on fibronectin-coated coverslips (BD Biosciences) at a density of 2 x 10° 
cells/cm?. After incubation at 37°C for 48 hours, cells were infected with IBV M41, and at the 
indicated time points fixed in 70% ethanol at 4 *C. Hybridisation of the probes was performed 
using the manufacturer's protocol for adherent cells. Imaging was performed using a Roper (Evry, 
France) Spinning Disc Confocal System on a Nikon Eclipse Ti microscope using a 100 x Plan apo oil 
immersion objective (NA 1.4) and a 491 nm laser line. Z-stacks were collected with 0.25 um Z- 


intervals. For each channel, maximum Z-stack projections were made and processed with ImageJ. 


Chicken type 1 IFN bioassay. 

Bioactive chicken type I interferon (chIFN) was measured using a bioassay based on the CEC-32 
quail reporter cell line expressing luciferase under the control of the chicken Mx promotor (47) 
(kindly provided by Prof. Peter Staeheli). Briefly, CEC-32 were incubated with serial dilutions of 
chIFN-containing samples for 6 hours, after which luciferase activity was quantified and IFN 
concentrations calculated using a chIFN standard. To avoid influence of IBV on the assay, samples 


were heat inactivated at 56°C for 30 min, which did not influence type I chIFN bioactivity. 


Gene silencing 

siRNAs targeting chicken 7/r3 and Mda5 were designed by and purchased from Microsynth, 
Switzerland (sequences in Table 1). Transfections were performed using siLentFect (Biorad) at a 
final siRNA concentration of 20 nM. For one well, 160 ng siRNA was combined with 1 uL siLentFect 
in 100 uL OptiMEM (Gibco) and incubated for 20 min. The siRNA complexes were added to 2 x 10" 
DF-1 cells grown in 500 uL medium per well in a 24 well plate. siRNA complexes were left on the 


cells for 48 hours before further experiments were performed. 
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Statistics 

All statistical analyses were performed in GraphPad Prism 5.0. RT-qPCR fold changes were first log 
transformed and then used for statistical analysis. For all tests, equality of variance was assessed 
using Bartlett's test. Significant differences (P<0.01) were determined by a one way or two-way 


ANOVA (indicated in the figure legend) followed by a Bonferroni post-hoc test. 





233 


234 


235 


236 


237 


238 


239 


240 


241 


242 


243 


244 


245 


246 


247 


248 


249 


250 


251 


252 


253 


254 


255 


256 


257 


258 


259 


260 


261 


262 


263 


264 


Results 


IBV delays the onset of an IFN response during infection of primary chicken cells 

To investigate the kinetics of viral replication and IFN induction upon infection with the avian 
Gammacoronavirus IBV, we infected primary CEK cells (24) with the IBV M41 strain. To monitor 
the kinetics of the IFN response in relation to IBV replication, we quantified transcription of IfnB, a 
set of genes involved in innate immunity, extracellular IFN protein production, virus titres and IBV 
RNA in M41-infected CEK cells. In line with previous observations (48), progeny virus was produced 
after 6 hpi and virus titres reached a maximum around 24 hpi (Fig. 1A). Total intracellular IBV RNA 
levels reflected the kinetics of infectious IBV virus in the supernatant (Fig. 1A), reaching maximum 
levels around 24 hpi. IfnB expression was delayed with respect to the peak of viral replication and 
remained low until 18 hpi, after which it was strongly upregulated, peaking around 36 hpi (Fig. 
1B). IFN protein activity levels were quantified using a chicken IFN-specific Mx-luc cell-based 
bioassay showing accumulation of IFN from 36 hpi onwards (Fig. 1B). Concomitant with IfnB, a 
subset of genes involved in innate immunity, including Mx, Oas and 1/8, were upregulated whereas 
others, such as Tir7, Adar, Isg12, MHC-I and Ifnar2 appeared not, or only marginally affected by 
IBV infection (Fig. 1C). Pattern recognition receptors Mda5 and TIr3 and the transcription factor 
Irf3 were also upregulated (Fig. 1C), which is of interest given the role of these PRRs in virus 


recognition. 


The delayed IFN response is independent from the cell type or virus strain 

Ifnf transcription during infection with coronaviruses such as MHV and SARS-Co is generally low 
(9, 10, 12, 13, 49), and was shown to be dependent on cell type and virus strain (50). The delayed 
induction of IfnB transcription observed in IBV M41-infected CEK cells prompted us to investigate 
whether induction of IfnB would be dependent on the cell type or IBV strain. Epithelial cells isolated 
from trachea of 10-week-old SPF chickens and DF-1 chicken fibroblast cells were infected with IBV 
M41 or IBV Beaudette (Beau-R, (38)). At several time points after infection, IfnB levels were 
monitored by RT-qPCR (Fig. 2A and 2B). In both cell types IfnB transcription followed the same 
kinetics observed in CEK cells (Fig. 1B), indicating that induction of IfnB by IBV is independent of 
cell type. To study whether induction of IfnB transcription differs between different strains of IBV, 
we also infected CEK cells with the QX and ItO2 strains of IBV (Fig. 2C). Although we observed 


some differences in absolute levels of IfnB upregulation induced by OK, ItO2 and M41, kinetics of 
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IfnB transcription were similar, suggesting that delayed induction of IfnB transcription could be 
considered a general feature of IBV infection in chicken cells. 

To assess whether CEK and DF-1 cells do have the intrinsic ability to express IfnB earlier than 18h, 
we stimulated these cells with extracellular polyI:C (pl:C), transfected pI:C (t[pI:C]) or with the 
dsRNA virus Infectious Pancreatic Necrosis Virus (IPNV). We found that stimulation of CEK cells 
with pI:C could induce IfnB transcription as early as one hour after stimulation (Fig. 3A). In DF-1 
cells, stimulation with IPNV, and t[pI:C], but not pI:C, induced IfnB already at 4h (Fig. 3B). The 
observation that DF-1 cells do not respond to stimulation with extracellular dsRNA, is in accordance 
with previous findings and is most likely due to the lack of surface expression of TLR3 (51). In 
addition, a 12h infection of CEK cells with Sindbis, IPNV or Rift Valley Fever Virus clone 13 (RVFV 
CI13) induced a clear transcription of IfnB (Fig. 3C). These results suggest that delayed expression 


of IfnB is a specific feature of IBV infection and not an intrinsic characteristic of chicken cells. 


The intracellular pattern recognition receptor MDAS is the primary sensor of IBV 

In general, dsRNA has been shown to be the canonical inducer of IfnB during infection with Alpha- 
and Betacoronaviruses (16, 50). To determine which pattern recognition receptor (PRR) would be 
involved in sensing (ds)RNA of the Gammacoronavirus IBV, leading to subsequent IfnB 
transcription, we first examined the possibility that IBV-(ds)RNA could be sensed extracellularly by, 
for example, cell-surface receptors. To investigate this, CEK cells were infected with IBV M41 in the 
presence of RNase A and IfnB expression was analysed. As a positive control, CEK cells were 
stimulated with pI:C in the presence or absence of RNase A. The IFN response to pI:C was greatly 
inhibited by addition of RNase A, which had no effect on Ifnf levels induced by infection with IBV 
M41 (Fig. 4A). These data suggest that IfnB upregulation during the late stage (>18 hpi) of IBV 
infection could be the result of sensing of IBV-(ds)RNA by an intracellular rather than an 
extracellular pattern recognition receptor. This is consistent with our observation that IBV infection 
can be detected by DF-1 cells, which show only a marginal upregulation of IfnB transcription in 
response to extracellular dsRNA (see Fig. 3B). In general, dsRNA can be recognised by membrane- 
bound TLR3 and cytosolic RLRs such as MDA5 and RIG-I. Genome mining strongly indicates that 
chickens do not express a RIG-I homologue (52), leaving TLR3 and MDA5 as the two PRRs 
potentially involved in dsRNA sensing. Silencing of Mda5, but not 7/r3, in DF-1 cells resulted in a 
70% decrease in IfnB transcription (Fig. 4B). Similar results were obtained with an Ifnf-luc DF-1 


reporter cell line in which silencing of Mda5, but not 7/r3, resulted in a 70% decrease in luciferase 
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activity by the reporter cells (Fig. 4C). Because no antibody against chicken MDAS is currently 
available for protein detection, successful knockdown was evaluated using RT-qPCR demonstrating 
a silencing efficiency for both T/r3 and Mda5 of approximately 60% (data not shown). Replication of 
IBV at the investigated time point was not affected by knockdown of neither 7/r3 nor Mda5, as 
measured by both virus titre and intracellular IBV total RNA (Fig. 4D). These results indicate that 


MDAS is the primary PRR responsible for sensing Gammacoronavirus IBV-(ds)RNA in chicken cells. 


Early accumulation of dsRNA in IBV-infected cells does not result in early induction of 
Ifnf 

Having assessed that chicken cells can indeed promptly respond to stimulation with dsRNA (Fig. 3) 
and having identified MDA5 as the primary sensor involved in the detection of IBV (Fig. 4), we 
investigated whether there would be a temporal difference between IBV-induced accumulation of 
dsRNA and the upregulation of IfnB transcription in CEK cells. Indeed, dsRNA could clearly be 
detected, even at low MOI of 0.01, by 12 hpi (Fig. 5A). In contrast, IfnB levels at this time point 
remained low (Fig. 5B) even in cell cultures infected at higher MOIs of 1 or 10 and despite the 
increased abundance of dsRNA. To further investigate the time lag between early accumulation of 
dsRNA and late IfnB expression, we performed a time course analysis. Foci of dsRNA could be 
detected as early as 3 hours post-infection only in IBV-infected cells (Fig. 5C, inset 3hpi) indicating 
that dsRNA starts accumulating very early in IBV-infected cells but apparently only leads to late 
(>18 hpi) IfnB transcription. 

Primary CEK cells consist of a heterogeneous mix of cell types. Even at high MOI, IBV M41 infects 
only “70% of the cells, indicating that not all cells are permissive to IBV M41 infection. In order to 
assess whether the time lag between accumulation of dsRNA and IfnB expression could be due to 
the induction of IfnB in bystander rather than IBV-infected cells, we used RNA fluorescent in situ 
hybridisation to visualise IfnB mRNA in IBV-infected CEK cell cultures (Fig. 5D). At 12 hpi and low 
MOI (0.1), with most cells showing clear foci of dsRNA, none of the IBV-infected cells displayed an 
accumulation of IfnB mRNA. At 12 hpi and a higher MOI, a few cells stained positive for IfnB mRNA 
and only later, at 24 hpi, did most IBV-infected cells also stain positive for IfnB mRNA, the kinetics 
of which closely following that observed in Figure 5A. In all cases, detection of IfnB mRNA was 
restricted to cells that contained dsRNA. Altogether our data shows that IBV-infected, but not 


adjacent uninfected cells, upregulate IfnB transcription in response to IBV infection. The significant 
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time lag between accumulation of dsRNA and Ifnf transcription further suggests the presence of a 


mechanism adopted by IBV to circumvent the onset of an IFN response. 


Accessory proteins 3a and 3b regulate IFN transcription and protein production 

To investigate whether the accessory proteins of IBV might play a role in the observed delay in IfnB 
transcription, we infected CEK cells with IBV scAUG3ab and scAUG5ab null viruses and the parental 
Beau-R virus (scAUG viruses possess a scrambled AUG start codon resulting in transcription but not 
translation of either ORFs 3a and 3b or 5a and 5b (36, 37)). Infection with the scAUG3ab, but not a 
scAUG5ab null virus resulted by 24 hpi in increased upregulation of IfnB expression (Fig. 6A). 
Indicating that either one, or a combination of, accessory proteins 3a and 3b play a role in down 
regulating IfnB transcription. The difference in IfnB transcription between the scAUG3ab and the 
parental (Beau-R) virus could not be ascribed to differences in kinetics of virus replication, as all 
viruses displayed similar growth kinetics until 24 hpi (Fig. 6B). To determine whether 3a, 3b or 
both accessory proteins are involved in the observed down regulation of the IFN response, we 
quantified IfnB transcription and IFN protein production in CEK cells infected with scAUG3a, 
SCAUG3b and scAUG3ab mutant viruses, and compared the values observed in cells infected with 
Beau-R (Fig. 6C and 6D). Infection with all mutant viruses led to an increased transcription of IfnB 
when compared to the Beau-R (Fig. 6C), indicating that the presence of either one of the two 
accessory proteins is sufficient to limit IfnB transcription. The kinetics of IfnB transcription in 
response to AUG3a/b differs between Fig. 6A and Fig. 6C. In Fig 6C there is a significant difference 
in IfnB transcription between AUG3a/b and Beau-R, which is absent in Fig 6A. This difference can 
probably be attributed to variation in the kinetics of IfnB transcription between primary CEK cells 
isolated from embryos originating from different flocks. Nonetheless, this difference does not affect 
the conclusion that knockout of 3a and 3b leads to an increase in transcription of IfnB. 
No significant differences in IFN protein production were observed between cells infected with the 
Beau-R and the scAUG3ab double null virus, except at 36 hpi. However, infection with scAUG3b 
virus led to an increase in IFN protein levels, whereas infection with the scAUG3a virus led to a 
decrease in IFN when compared to both Beau-R and scAUG3ab double null virus (Fig. 6D). Taken 
together, these results indicate that accessory proteins 3a and 3b both play a role in the inhibition 
of IfnB transcription but have distinct and opposing effects on protein production. Accessory protein 


3b seems to be involved in limiting IFN protein activity whereas 3a is involved in promoting it. 
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Signalling of non-self dsRNA remains intact in IBV-infected cells 

Since IBV showed the intrinsic ability to delay IfnB transcription in several cell types (Fig. 1 and 
Fig. 2) even in the presence of high levels of intracellular dsRNA (Fig. 5), we investigated the 
ability of IBV to interfere with sensing of non-self dsRNA by TLR3 or MDA5. We infected CEK cells 
with IBV M41 and subseguently used extracellular poly I:C to trigger TLR3 signalling (Fig. 7A). 
Stimulation with pI:C alone led to a significant increase in IfnB transcription whereas stimulation 
with pI:C following an infection with IBV led to an enhanced increase in IfnB transcription in an 
MOI-dependent manner. These results indicated that IBV infection does not interfere with TLR3- 
mediated IfnB transcription, on the contrary IBV infection appears to result in a synergistic 
activation of the TLR3 pathway triggered by pI:C. Next, we investigated whether IBV infection 
could interfere with MDA5-mediated transcription of IfnB. Although transfection of pI:C into the 
intracellular compartment is a commonly used ligand of MDA5, this method induced very little 
transcription of IfnB in primary CEK cells, because of low transfection efficiency (data not shown). 
As an alternative route to stimulate MDA5 in primary chicken cells, we investigated the use of 
either RVFV CI13 or IPNV, that induce IfnB transcription in CEK cells (Fig 3C). RVFV Cl13 is a (-) 
ssRNA virus with a truncated IFN antagonist (53), for which RIG-I, but not by MDA5 or TLR3, was 
previously shown to be the most likely PRR in mammalian cells (54, 55). Since chickens, as 
opposed to most mammals, do not have a RIG-I homologue, the most likely PRR for RVFV in CEK 
cells would be MDAS. IPNV is a birnavirus with a dsRNA genome that naturally infects salmonids 
but has been shown to enter but not replicate in cells of warm-blooded animals (56). To date, the 
PRR responsible for sensing IPNV dsRNA has not been described. Knockdown experiments in DF-1 
IfnB-luc reporter cells, using siRNAs against chicken MDA5 or TLR3, revealed that MDA5, but not 
TLR3, is the prime PRR for IPNV (Fig. 7B). These findings were confirmed using MEFs (mouse 
embryo fibroblasts) from knockout mice deficient in expression of either MDA5, RIG-I or the 
downstream adaptor protein MAVS. Here, knockout of either MDA5 or MAVS abrogated sensing of 
IPNV as shown by a strong reduction of IfnB transcription, whereas knockout of RIG-I did not (Fig. 
ፖር). Both IPNV and RVFV Cl13 were subsequently used to investigate whether IBV infection could 
interfere with MDA5-mediated transcription of Ifnf in CEK cells. 

Using quantification of IfnB transcription by RT-qPCR as read out, we could show that IBV infection 
does not interfere with MDA5-mediated signalling of IPNV (Fig. 7D) or RVFV Cl13 (Fig. 7E), in fact 
it had a synergistic effect on IfnB transcription as previously observed for TLR3-mediated signalling 


(Fig. 7A). Similar results were obtained when stimulating IBV-infected DF-1 cells with IPNV or 
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t[pI:C] (Fig. 7F), indicating that the observed synergistic effect is not specific to CEK cells. Taken 
together, IBV infection very efficiently prevents sensing of IBV (ds)RNA, but our results indicate 


that it does not interfere with sensing and downstream signalling of other non-self (ds)RNA ligands. 


DISCUSSION 

In this study we performed a comprehensive analysis of the kinetics of IBV infection in avian cells 
and studied the mechanisms by which IBV interferes with the onset of the type ፤ IFN response. We 
show that infection with the Gammacoronavirus IBV leads to a considerable activation of the type I 
IFN response, albeit delayed with respect to the peak of viral replication and accumulation of viral 
dsRNA. Using an siRNA knockdown approach we show that MDA5 is the main receptor involved in 
the induction of IfnB expression during IBV infection. We present evidence that IBV accessory 
proteins 3a and 3b play a role in the modulation of the delayed IFN response, by regulating 
interferon production both at the transcriptional as well as translational level. In addition, we show 
that although IBV alone effectively prevents IfnB induction in IBV-infected cells, it does not block 
IfnB induction upon stimulation of IBV-infected cells with other RIG-I, MDA5 or TLR3 ligands. To 
our knowledge, this study provides the most comprehensive analysis of the interplay between a 


Gammacoronavirus and the avian type 1 IFN response. 


Much of our current knowledge about the interaction of coronaviruses with the innate immune 
response (reviewed in (57)) comes from studies in mice and mouse cells using mouse hepatitis 
virus (MHV). MHV activated IFN production only in specific cell types and an efficient IFN response 
was only mounted by plasmacytoid dendritic cells (58), bone marrow derived macrophages (10, 
59) and oligodendrocytes (10). In a recent study on SARS-CoV and MERS-CoV in an epithelial lung 
cell line, ISGs started to be upregulated at 12 hpi (60), when virus titres were already reaching 
their maximum. The kinetics of IFN response observed in our study are in line with aforementioned 
studies, however it must be noted that in most cell types, infection with A/pha or Betacoronaviruses 
induced very little, if any, IfnB transcription (8-13, 49). This suggests that all coronaviruses are 


able to modulate the activation of the type I IFN response. 


We found that IBV infection is detected by various chicken cell types, but until now it was unknown 
which PRR was involved. MHV has been shown to be detected by MDA5 and not RIG-I or TLR3 in 


brain macrophages (50), by both MDA5 and RIG-I in an oligodendrocyte derived cell line (10) and 
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by TLR7 in plasmacytoid dendritic cells (58). Analysis of the chicken genome suggests that chicken 
lack a RIG-I homologue (52), and basal expression of T/r7 was found to be very low in CEK cells 
(data not shown). We therefore silenced the remaining candidate RNA sensors MDA5 and TLR3, 
and were able to show that MDA5, but not TLR3, is involved in the sensing of IBV. Silencing of 
Mda5 did not lead to an increase in replication of IBV, suggesting that IBV might have developed 


strategies to counteract the activated IFN response. 


We recently reported membrane rearrangements in chicken cells infected with IBV (48), similar to 
those found in cells infected with Betacoronaviruses. In theory, the formation of intracellular 
membrane rearrangements might partly explain the discrepancy observed in the kinetics of dsRNA 
accumulation and IfnB upregulation. Indeed, for SARS-CoV it has been shown that virus-induced 
double membrane vesicles (DMVs) contain dsRNA (14), suggesting that coronaviruses might 
exploit membrane structures to shield dsRNA from recognition by host PRRs (61). However, the 
kinetics of IfnB transcription were not investigated in these studies. The presence of coronavirus- 
induced DMVs has been demonstrated as early as 2 hpi in SARS-CoV-infected cells (14). Although 
we did not demonstrate the presence of DMVs in IBV-infected chicken cells at time points earlier 
than 7 hpi (48), it is likely that DMVs could also be present at earlier time points. As such, the 
timing of DMV formation in coronavirus-infected cells could suggest that membrane 
rearrangements play a role in the delayed activation of the IFN response by shielding dsRNA from 


cellular PRRs. 


In addition to membrane rearrangements, coronavirus-encoded proteins, including numerous 
accessory genes, have been shown to interfere with the type I IFN response pathway (reviewed in 
(17, 18). To investigate the possible role of IBV accessory proteins in the regulation of the IFN 
response we made use of our previously constructed mutant IBV Beau-R viruses that do not 
express either one or more of the four accessory proteins 3a, 3b, 5a and 5b. Previously we have 
demonstrated the accessory genes of IBV are not essential for replication (36, 37). In the present 
study we show that infection of CEK cells with 3a or 3b null viruses as well as a 3a/3b double null 
virus led to increased IfnB transcription compared to Beau-R. Because the kinetics of IfnB 
transcription of 3a, 3b and 3a/3b null viruses are comparable to the parental virus, we conclude 
that, 3a and 3b are probably not responsible for the delay in IfnB transcription, suggesting that IBV 


utilises additional strategies to delay transcription of IfnB. Apart from their effect on IfnB 
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transcription, 3a and 3b seem to have opposing effects on IFN protein production by IBV infected 
cells. Infection with the 3b null virus resulted in increased IFN production whereas infection with 
the 3a null virus resulted in reduced IFN levels compared to the Beau-R virus. Together with the 
observation that IFN production induced by the 3a/3b double null virus is comparable to that 
induced by Beau-R virus, our data suggests that accessory proteins 38 and 3b antagonise each 


other to tightly regulate IFN production (Fig. 6B). 


Using the eukaryotic linear motif server (62), we identified a Protein phosphatase 1 (PP1)-binding 
UKISF? domain in the IBV 3b protein sequence. The canonical PPi-binding motif is 
[R/K][V/I/L]X[F/W], in which x can be any amino acid except proline (63). Interestingly, 
Alphacoronavirus TGEV accessory protein 7 (TGEV-7) has been shown to bind PP1 via a binding 
motif similar to that found in IBV 3b (64). Similar to IBV scAUG3b, infection with TGEV-A7 led to 
increased mRNA and protein levels of IFNB (65). The fact that both TGEV-7 and IBV 3b contain a 
PP1 binding domain indicates that interaction with PP1 could be a common strategy of 
coronaviruses to inhibit the host innate immune response. The mechanism by which interaction of 
coronavirus accessory proteins with PP1 counteracts the innate immune response still needs to be 
determined. One clue might come from the PPi-binding domain of Measles virus V, which was 
recently shown to be essential for inhibition of MDAS signalling (66, 67). Measles V protein binds 
PP1 and inhibits dephosphorylation of MDA5, which is required for activation and subsequent 
signalling by MDA5. Motif analysis for IBV 3a protein did not reveal the presence of relevant 
motives that might explain the observed activity of 3a on IFN regulation. We conclude that both 
accessory proteins 3a and 3b limit IfnB transcription but have distinct and opposing effects on 
protein production. Whilst 3a seems to promote IFN production, 3b seems to be involved in limiting 
IFN protein production, possibly through a similar mechanism as described for protein 7 of TGEV. 
The fact that IBV 3a and 3b have opposing roles in regulating IFN production indicates that CoV's 
tightly regulate IFN production to balance their own survival with that of the host. This hypothesis 
is supported by the observation that field isolates lacking 3a and 3b display reduced virulence in 
vitro as well as in vivo (68). Elucidation of the exact mechanisms of action of 3a and 3b will be the 


subject of further investigation. 
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To investigate whether IBV interferes with a general sensing of (ds)RNA ligands or downstream 
signalling that leads to IfnB transcription, we stimulated IBV-infected cells with TLR3, RIG-I, and 
MDAS ligands. Surprisingly, we found that infection with IBV did not reduce IfnB transcription but 
rather increased IfnB levels upon stimulation with these PRR ligands. Similar to IBV, MHV has been 
shown unable to inhibit expression of IfnB induced by either t[pI:C] or Sendai virus (69, 70), but in 
these studies no synergistic effect was observed. Currently, we can only speculate about the cause 
of this synergistic effect. It appears that IBV infection ‘arms’ the IfnB induction pathway, without 
actually triggering it, possibly by enhancing the activity of one or more components of the pathway 
leading to IfnB upregulation. One possibility is that IBV-proteins interact with host-proteins that 
regulate this pathway through ubiquitination and phosphorylation (reviewed in (71)). The fact that 
stimulation with either TLR3 or MDAS ligands resulted in exacerbated transcription of IfnB indicates 


that IBV influences a component which is downstream of both MDA5 and TLR3. 


Taken together, our study provides the first comprehensive analysis of host-virus interactions of a 
Gammacoronavirus with the avian innate immune response. We show that the Gammacoronavirus 
IBV, induces activation of the type I IFN response in primary chicken renal cells, tracheal epithelial 
cells and in a chicken cell line. We show that activation of the IFN response is dependent on MDA5 
but is delayed with respect to the peak of virus replication. We demonstrate that IfnB transcription 
is restricted to IBV-infected, dsRNA-containing cells and provide evidence that accessory proteins 


3a and 3b of IBV are involved in regulating transcription as well as protein production of type I IFN. 
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FIGURE LEGENDS 


FIG 1. IBV infection delays Ifnf upregulation. 

Chicken embryo kidney (CEK) cells were infected with IBV M41 at an MOI of 0.1. (A) Replication of 
IBV was quantified by titration in cell culture supernatants of infected cells; in a parallel 
experiment, intracellular IBV RNA was quantified using RT-gPCR. (8) IfnB mRNA levels were 
determined using RT-gPCR and IFN protein levels using a chicken IFN-specific mx-luc cell-based 
bioassay, respectively. (C) Expression of genes involved in the antiviral response. All gene 
expressions were calculated as fold changes relative to uninfected control cells and normalised 
against an external reference gene (luciferase). For IBV total RNA, fold changes were calculated 
relative to C, 30. Depicted are the results of a representative experiment out of three independent 


experiments. 


FIG 2. Delayed induction of Ifnf transcription is independent of cell type or IBV strain. 

(A) Epithelial cells from adult chicken trachea were infected with IBV M41 at MOI 0.1. (B) 
Fibroblast DF-1 cells were infected with IBV Beau-R at MOI 0.1. (ርን CEK cells were infected with 
IBV strains M41, OK and It02, at an MOI 0.1. Intracellular IBV total RNA (open diamonds) and IfnB 
mRNA (bars) are depicted as fold changes as assessed by RT-gPCR. Gene expression of IfnB was 
calculated as fold changes relative to uninfected control cells and normalised against an external 
reference gene (luciferase).. For IBV total RNA, fold changes were calculated relative to C. 30. Error 


bars indicate standard deviation. 


FIG 3. Chicken cells have the intrinsic ability to respond rapidly to dsRNA. 

(A) CEK cells were seeded in 24 well plates and 48 hours later stimulated with extracellular poly 
I:C for the indicated times. (B) DF-1 cells were infected with IPNV, a non-replicating dsRNA virus, 
or stimulated with extracellular pI:C or transfected pl:C (t[pI:C]). Four hours later, IfnB fold 
changes were determined by RT-gPCR. Bars represent the mean (plus standard deviation) of 
triplicate wells from a representative experiment. Asterisks indicate significant differences (P<0.01) 
with respect to the non-stimulated control as assessed by one-way ANOVA followed by a Bonferroni 
post-hoc test. (ርን CEK cells were infected with IBV M41 (MOI 1), IBV Beau-R, (MOI 1), Sindbis-GFP 
(MOI 1), IPNV (MOI 50) and RVFV Cl13 (MOI 5). Depicted are IfnB fold changes at 12 hpi relative 


to uninfected control cells as assessed by RT-gPCR. 
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FIG 4. MDA5, and not TLR3, is the prime sensor of IBV. 

(A) CEK cells were infected with IBV M41 for 24 hours, in the presence or absence of RNase A. IfnB 
expression was analysed by RT-gPCR. Stimulation with pI:C in the presence or absence of RNase A 
was included as a positive control. (B,D) DF-1 cells and (C) DF-1 IfnB-luc reporter cells were 
transfected with siRNAs against 7/r3, Mda5 or a control siRNA and 48 hours later infected with IBV 
M41 (MOI 0.1). (B) IfnB mRNA, (ር) IfnB-luciferase activity, and (D) IBV titres and intracellular RNA 
were analysed 18 hpi. Bars represent the mean (plus standard deviation) of triplicate wells from a 
representative experiment. Asterisks indicate significant differences (P<0.01) with respect to the 
non-RNaseA-treated control (A) or to the siRNA control (B-C), as assessed by one-way ANOVA 


followed by a Bonferroni post-hoc test. 


FIG 5. Early accumulation of dsRNA in IBV-infected cells does not result in early 
induction of IfnB 

CEK cells were infected with IBV M41 or IBV Beau-R at the indicated MOIs. At time point 12 hpi (A) 
dsRNA was visualised in M41-infected cells using an antibody against dsRNA. (B) Expression of IfnB 
mRNA was analysed by RT-gPCR. (ር) CEK cells were infected with IBV M41 and accumulation of 
dsRNA was visualised at the indicated time post infection. (D) RNA fluorescent in situ hybridisation 
of IfnB mRNA in IBV M41-infected CEK cells. Open arrowheads indicate cells that contain dsRNA 


and no IfnB mRNA. Solid white arrowheads indicate cells that contain both dsRNA and IfnB mRNA. 


FIG 6. Accessory proteins 3a and 3b are involved in regulation of IFN transcription and 
protein production. 

(A) CEK cells were infected with IBV Beau-R 3a/3b (scAUG3ab) or 5a/5b (scAUG5ab) null viruses 
(MOI 0.1). IfnB levels were determined using RT-gPCR. (B-D) CEK cells were infected with 
scAUG3a, scAUG3b or scAUG3ab null IBV viruses (MOI 0.1). In the same cultures (B) IfnB mRNA, 
(ር) virus titres and (D) type I IFN protein were guantified. Bars represent the mean (plus standard 
deviation) of triplicate wells from a representative experiment. Significant differences (P<0.01) 
relative to the Beau-R virus at the same timepoint (ችን or between the indicated bars (#) as 


assessed by a two-way ANOVA followed by a Bonferroni post-hoc test. 
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FIG 7. Signalling of non-self RNA remains intact in IBV-infected cells 

(A) CEK cells were infected with IBV M41 for 3 hours and stimulated with extracellular poly 1:C (50 
ug/ml) for an additional 3 hours after which Ifnf transcription was analysed by RT-gPCR. (B) DF-1 
Ifnf-luc reporter cells were transfected with siRNAs against 7/r3, Mda5 or a control siRNA and 48 
hours later infected with IPNV (MOI 50); at 6 hpi luciferase activity was quantified. (C) Knockout 
(KO) and wild-type (wt) MEFs were infected with IPNV (MOI 50) for 8 hours. (D) CEK cells were 
infected with IBV M41 (MOI 10) for 6h and super-infected with IPNV or UV-inactivated IPNV (MOI 
50) for an additional 6 h. (E) CEK cells were co-infected with IBV M41 (MOI 5) and RVFV clone 13 
(MOI 5) and sampled at 6 hpi. (F) DF-1 cells were infected with IBV Beau-R (MOI 1) for 3 h and 
super-infected with IPNV (MOI 50) or transfected with pI:C (t[pI:C], 500 ng/well) for an additional 
4 h. (C-F) IfnB levels were quantified by RT-qPCR. Bars represent the mean (plus standard 
deviation) of triplicate wells. Significant differences (P«0.01) are indicated by (*) as assessed by 


one-way ANOVA followed by a Bonferroni post-hoc test. 
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TLR3 siRNA 5 UCGAAUACUUGGCUUUAAA NM_001011691 
AS UUUAAAGCCAAGUAUUCGA 

ctrl siRNA 5 AGGUAGUGUAAUCGCCUUG - 
ለ5 CAAGGCGAUUACACUACCU 

MDA5 siRNA 5 ACACUGGUAUCAAGUUAUU GU570144 
AS AAUAACUUGAUACCAGUGU 

IFNB RO primer FW  GCTCTCACCACCACCTTCTC ENSGALTO0000039477 
RV GCTTGCTTCTTGTCCTTGCT 

IFNa RQ primer FW ATCCTGCTGCTCACGCTCCTTCT XM_004937096 40 
RV GGTGTTGCTGGTGTCCAGGATG 

IRF3 RQ primer FW CAGTGCTTCTCCAGCACAAA NM_205372 
RV  TGCATGTGGTATTGCTCGAT 

IRF1 RO primer FW CAGGAAGTGGAGGTGGAGAA ENSGALG00000006785 
RV  TGGTAGATGTCGTTGGTGCT 

TLR7 RQ primer FW TTCTGGCCACAGATGTGACC NM 001011688 40 
RV CCTTCAACTTGGCAGTGCAG 

TLR3 RO primer FW TCAGTACATTTGTAACACCCCGCC NM 001011691 40 
RV GGCGTCATAATCAAACACTCC 

MDA5 RO primer FW TGGAGCTGGGCATCTTTCAG 6ህ570144 
RV GTTCCCACGACTCTCAATAACAGT 

Mx RQ primer FW TTGTCTGGTGTTGCTCTTCCT ENSGALTO0000025999 
RV GCTGTATTTCTGTGTTGCGGTA 

OAS RQ primer FW CACGGCCTCTTCTACGACA NM_205041 41 
RV TGGGCCATACGGTGTAGACT 

IL8 RQ primer FW TTGGAAGCCACTTCAGTCAGAC NM_205498 41 
RV GGAGCAGGAGGAATTACCAGTT 

PKR RQ primer FW CCTCTGCTGGCCTTACTGTCA NM_204487 42 
RV  AAGAGAGGCAGAAGGAATAATTTGCC 

ADAR RQ primer FW TGTTTGTGATGGCTGTTGAG AF403114 
RV  AGATGTGAAGTCCGTGTTG 

[5612 RO primer FW TAAGGGATGGATGGCGAAG NM 001002856 
RV GCAGTATCTTTATTGTTCTCAC 

MHC-I RO primer FW  CTTCATTGCCTTCGACAAAG NM_001031338 41 
RV GCCACTCCACGCAGGT 

IFNAR2 RQprimer FW  GCTTGTGTICGTCAGCATT ENSGALT00000036778 41 
RV  TTCGCAATCTTCCAGTTGT 

IBV-N RQ primer FW GAAGAAAACCAGTCCCAGA AY851295 
RV TTACCAGCAACCCACAC 

Luciferase RO primer FW TGTTGGGCGCGTTATITATC X65316 


RV AGGCTGCGAAATGTTCATACT 
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